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Dilcpton and photon production in heavy-ion collisions at SPS energies is studied 
in the relativistic transport model that incorporates self-consistently the change of 
hadron masses in dense matter. It is found that the dilepton spectra in proton- 
nucleus reactions can be well described by the conventional mechanisms of Dalitz 
decay, primary vector meson decay, decay of charmed mesons, and the initial 
Drell-Yan processes. However, to provide a quantitative explanation of the ob- 
served dilepton spectra in central heavy-ion collisions requires contributions other 
than these direct decays and also various medium effects. Introducing a decrease 
of vector meson masses in hot dense medium, we find that the low-mass dilep- 
ton enhancement can be satisfactorily explained. Furthermore, to explain the 
intermediate-mass dilepton enhancement in heavy-ion collisions, secondary pro- 
cesses such as nai — > 11 are found to be very important. Finally, the single photon 
spectra in our calculations with either free or in-medium meson masses do not 
exceed the upper limit measured by the WA80 Collaboration. 



1 INTRODUCTION 

The experimental measurement and theoretical investigation of electromag- 
netic observables in heavypion collisions constitutes one of the most active and 
exciting fields in physics u. Because of their relatively weak final-state in- 
teractions with the hadronic environment, dileptons, as well as photons, are 
considered ideal probes of the early stage of heavy-ion collisions, where quark- 
m-plasma (QGP) formation and chiral symmetry restoration are expected 
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Dilcpton mass spectra can be basically divided into three mass regions. 
The low-mass region below to^ (~ 1 GeV) is dominated by hadronic interaction 
and hadronic decay in the freeze-out. In the intermediate-mass region between 
m<£ and about 2.5 GeV, the contribution from the thermalized QGP might be 
seen. In the high-mass region at and above mj/^ the major effort has been 
the detection and understanding of J/\& suppression. 

So far, the experimental measurement of dilepton spectra in heavy-ion col- 
lisions has mainly been carried out at the CERN SPS by three collaborations: 
the GERES collaboration is dedicated to dielectron spectra in the low-mass re- 
gionEffi, the HELIOS-30 collaboration has measured dimupn spectra from the 
its threshold up to the J/ty region, and the NA38/NA50EI collaboration mea- 
sures dimuon spectra in the intermediate- and high-mass region, emphasizing 
the J/ 1 ^ supression. 

Recent observation of the enhancement of low-mass .dileptons in central 
heavy-ion collisions by the CERES clta and the HELIOS-30 collaborations has 
generated a great deal of interest in the heavy-ion community. Different dy- 
namical models have been used to calculate and understand these data. The 
results from many groups with standard scenarios (i.e., using vacuum meson 
properties) are in remarkable agreement with each other, but in significant 
disagreement with the data: the experimental^pecixa .in-the rnqgg region from 
fl.3-0.6 GeV are substantially vmderestimatedBBBByyfiiy (see also Ref. 
a). This has led to the suggestion of various medium effects that might be 
responsible for the obaaapd enhancement. In particular, the dropping vector 
meson mass scenario BOta due to the chiral symmetry restoration is found 
to provide a nice description of both the CERES and HELIOS-3 data. On 
the other hand, conventional many-body effects leading to a brodening of the 
rho-meson spectrum function also provide a good description of these datalia 

In the high-mass region around mjiy, the J/^Sl suppression has been a 
subject of great interest, since it was first proposed as a signal of the decon- 
fincmcnt phase transition to. Various investigations show that up to central 
S+Au collisions, the normal pre-resonance absorption in nuclear matter is suf- 
ficient to account for the observed J/^i> suppression. However, recent data 
from the NA50 collaboration for central Pb+Pb collisions show an additional 
strong 'anomalnus' suppression which might indicate the onset of the color 
deconfinement E3. 

Enhancement of dileptons in the intermediate-mass region from about 1 
GeV to about 2.5 GeV has also been observed by both the HELIOS-3 and 
NA38/NA50 collaborations in central S+W and S+U collisions as compared to 
that in the proton-induced reactions (normalized by the charged-particle multi- 
plicity) 0i3. Preliminary data from the NA50 collaboration also show significant 



2 



enhancement in central Pb+Pb collisions cl (see also Ref. □). The intermediate- 
mass dilepton spectra in heavy-ion collisions are particularly useful for the 
search of the QGP. It was originally suggested that in the intermediate-mass 
or the intermediate pt region, the electromagnetic xadiation from the QGP 
phase might shine over that from the hadronic phaseLirCJ. However, to extract 
from the measured dilepton spectra any information about the phase transi- 
tion and the properties of the QGP, it is essential that the contributions from 
the hadronic phase be precisely understood and carefully subtracted. 

Another piece of experimental data from CERN SPS that has been dis- 
posed extensively is the single photon spectra from the WA80 collaboration 
E3. In this experiment, tt° and r\ spectra were measured simultaneously, so 
that their contributions to the single photon spectra can be subtracted. It was 
found that the direct photon, excess over those background sources is about 5% 
for central S+Au collisionsE3. Similar measurement of inclusive photon spectra 
has been carried out by the CERES collaborationcil. The results are very much 
in agreement with those of WA80, namely, the inclusive photon spectra can 
be basically explained by hadronic decay at the freeze-out, particular radiative 
decay of 7r° and r\. 

In various hydrodynamics calculations the absence of significant thermal 



n of quark gluon plasma for- 
However, there is a major 



photon production has been used as an indk: 
mation that lowers the initial temperature £ 
uncertainty in this type of analysis, namely, the initial temperature depends 
sensitively on the degrees of freedom included in the analysis. Indeed, in Refs. 
I 13 H 27 1 it was found that including a sufficient amount of hadron resonances the 
WA80 photon data can be explained without-invoking the formation of the 
Detailed transport model calculations Ea support the findings of Refs. 



In this contribution, we will discuss (1) the enhancement of low-mass dilep- 
tons (Section 2), (2) the enhancement of intermediate-mass dileptons (Section 
3), and (3) the lack of direct photon signal (Section 4). A brief summary and 
outlook is presented in Section 5. 



2 Low-mass dilepton enhancement: a case for dropping rho meson 
mass 

To calculate dilepton spectra in heavy-ion collisions and to studying possible 
mediunx effects, the relativistic transport modelo based on the Walecka-type 
modelE3 has been quite useful, as it provides a thermodynamically consistent 
description of the medium effects through the scalar and vector fields. In heavy- 
ion collisions at SPS energies, many hadrons are produced in the initial nucleon- 
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nucleon interactions. This is usually modeled by the fragmentation of strings, 
which are the chromoelectric flux-tubes excited from the interacting quarks. 
One successful model—for taking into account this nonequilibrium dynamics 
is the RQMD modelEH. To extend the relativistic transport model to heavy- 
ion collisions at these energies, we have used as initial conditions the hadron 
abundance and distributions obtained from the string fragmentation in RQMD. 
Further interactions and decays of these hadrons are then taken into account 
as in usual relativistic transport model. 

To study the effects of dropping vector meson masses 00 on the dilepton 
spectrum in heavy-ion collisions, we have extended the Walecka model from 
the coupling of nucleons to scalar and vector fields to the coupling of light 
(quarks to these fields, using the ideas of the meson-quark coupling model 
c3. For a system of nucleons, pseudoscalar mesons, vector mesons, and axial- 
vector mesons at temperature T and baryon density ps, the scalar field (<j) is 
determined self-consistently from 



(2tt) 3 J E* N 



m„((T 
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eM(E* N - Hb)/T) + 1 exp((E* N + p B )/T) + 1 
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, dk^ — + / f, / dk—^- , . (1 ) 

2tt)3 J E* exp(E*/T) - 1 (2tt)3 J E* t exp{E*jT) -H> 

where we have used the constituent quark model relations for the nucleon and 
vector meson masses, i.e., m* N — tojv — 9<t(c)> m *p,u ~ m p,u — (2/3)g CT (<7), 
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the quark structure of the 77 meson in free space which leads to m* 
0A5g a (a), and the Weinberg sum rule relation between the rho-meson and a\ 
meson masses, i.e, m* ~ m ai la). 

The main contributions to dileptons with mass below 1.2 GeV are the 
Dalitz decay of 7r°, i] and u>, the direct leptonic decay of the primary p°, u> 
and </>, the pion-pion annihilation which proceeds through the p° meson, and 
the kaon-antikaon annihilation that proceeds through the <f> meson. The last 
two processes are unique in heavy-ion collisions, and especially the pion-pion 
annihilation is found to be very important for low-mass dilepton production. 
The explicit treatment of rho and phi meson formation, propagation, and decay 
in pion-pion and kaon-antikaon annihilation avoids possible double counting. 

The differential widths for the Dalitz decay of 7r°, 77, and oj are related to 
their radiative decay widths via the vector dominance model, which are taken 
from Ref. 123. The decay of a vector meson into the dilepton is determined by 
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the width, 



IW,- (M) = C7 I+I _ ^ (1 -_)/(! + _). (2) 

The coefficient C; +i - in the dielectron channel is 8.814 x 1CT 6 , 0.767 x 1(T 6 , 
and 1.344 x 10 -6 for p, u>, and <fi, respectively, and is determined from the 
measured width. For the dimuon channel, these values are slightly larger. The 
two-body processes (pion-pion and kaon-antikaon annihilation) are treated in 
two steps: the formation of the intermediate vector meson and its decay into 
dilcptons. 




Figure 1: Dilepton invariant mass spectra in p+Be (left window) and p+Au (right window) 
collisions at 450 GeV. 

The results for dilepton spectra from p+Be and p+Au collisions at 450 
GeV are shown in Fig. 1, together with data from the CERESH. It is seen that 
the data can be well reproduced by Dalitz decay of 7r°, r\ and u> mesons, and 
direct leptonic decay of p , to and tfi mesons. These results are thus similar to 
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those found in Rcf. El using the Hadron-String Dynamics and those constructed 
by the CERES collaboration from known and expected sources of dileptons 0. 
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Figure 2: Dilepton invariant mass spectra in S+Au collision at 200 AGeV (left window) and 
Pb+Au collisions at 160 AGeV (right window). The solid and dashed curves are obtained 
with in-medium and free meson masses, respectively. 

Our results for dilepton spectra in central S+Au collisions are shown in the 
left window of Fig. 2. The dashed curve is obtained with free meson masses. 
Although pion-pion annihilation is important for dilcptons with invariant mass 
from 0.3 to 0.65 GeV, it still does not give enough number of dileptons in this 
mass region. Furthermore, for masses around m PjW there are more dileptons 
predicted by the theoretical calculations than shown in_the experimental data. 
These are very similar to the results of Cassint el al. El based on the Hadron- 
String Dynamics model and Srivastava et al. El based on the hydrodynamical 
model. The results obtained with in-medium meson masses are shown by 
the solid curve. Compared with the results obtained with free meson masses, 
there is about a factor of 2-3 enhancement of the dilepton yield in the mass 
region from 0.2 to 0.6 GeV, which thus leads to a good agreement between the 
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theoretical results and the CERES data. Similar conclusions that dropping 
vector mespruxiasses can explain the CERES dilepton data have been obtained 
in Refs. 

In the right window of Fig. 2 we compare our prediction for central Au+Pb 
collisions at 160 AGeV with the preliminary data from the CERES collabo- 
ration. The normalization factor dN c h/drj here is the average charge particle 
pseudo-rapidity density in the pseudo-rapidity range of 2 to 3, and is about 
440 in this collision. In the results with free meson masses, shown by the 
dashed curve, there is a strong peak around m PjW , which is dominated by p 
meson decay as a result of an enhanced contribution from pion-pion annihi- 
lation in Pb+Au collisions than in S+Au and proton-nucleus collisions. In 
the case of in-medium meson masses, shown by the solid curve, the p meson 
peak shifts to a lower mass, and the peak around m p , w becomes a shoulder 
arising mainly from lo meson decay. At the same time we see an enhancement 
of low-mass dilcptons in the region of 0.25-0.6 GeV as in S+Au collisions. The 
agreement with the data is thus significantly improved when dropping vector 
meson masses are used. 

The same model has been used to calculate the dimuon spectra from p+W 
and central S+W collisions. The results for p+W are shown in the left window 
of Fig. 3. Again, the conventional mechanism describes quite well the measured 
spectra. The results for S+W collisions obtained with free meson masses are 
shown in the right window of Fig. 3 by the dashed curve, and are below 
the HELIOS-3 data in the mass region from 0.35 to 0.6 GeV, and slightly 
above the data around m PjW as in the CERES case. However, the discrepancy 
between the theory and the data is somewhat smaller in this case due to the 
smaller charged-particle multiplicity at a larger rapidity than in the CERES 
experiment. Our results obtained with in-medium meson masses are shown 
in the right window of Fig. 3 by the solid curve, and are in good agreement 
with the data. The importance of dropping rho meson_mass in explaining the 
HELIOS-3 data has also been found by Cassing et al. lij. 

It should be pointed out that, the effects of dropping rho meson mass 
on the low-mass dilepton spectra, as seen in Figs. 2 and 3, aiere predicted 
theoretically before the CERES data were published. In Ref. E-3, the RQMD 
model was used to calculate dilepton spectra in Au+Au collisions at AGS 
energies. It was found that, with an in-medium rho meson mass as predicted 
by the QCD sum rule calculation, the rho meson peak in the dilepton spectra 
shifts to about 0.5 GeV, and a significant enhancement is seen in this mass 
region. Similar observation was obtained in Ref. for dilepton spectra in 
heavy-ion collisions at SIS energies. 

Rho meson decays dominantly into two pions. Its spectral function can 
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Figure 3: Dilepton invariant mass spectra in p+W (left window) and S+W (right window) 
collisions at 200 AGeV. In the right window the solid and dashed curves are obtained with 
in-medium and free meson masses, respectively. 



therefore be modified by the medium modification of the pion dispersion func- 
tion. Furthermore, rho meson itself can interact with nucleons to form nucleon 
resonances, in particular iV(1720) and A(1905), which have appreciable decay 
branching ratios into the pN final state. Both effects lead to the broadening of 
the rho meson spectral function and thus the enhancement of low-mass dilert 
tons from the in-medium rho meson decay. This idea was shown in Ref. EU 
to also lead to a reasonable description of the low-mass dilepton enhancement 
observed by the CERES collaboration. 

3 Intermediate-mass dileptons: the importance of secondary pro- 
cesses 

To extend the previous calculation to the intermediate mass region, we need 
to include additional secondary processes in addition to irir and KK anni- 
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hilation, as well as the contributions from the charmed mesons and initial 
Drcll-Yan processes. The last two contributions involve hard processes which 
scale almost linearly with the participant nucleon number, and can thus be ex- 
trapolated from the proton-proton and proton-nucleus collisions. Such a study 
has recently been carried out by Braun-Munzinger et alE3. The results for the 
central S+W collisions corresponding to the HELIOS-3 acceptance are shown 
in the left window of Fig. 4, which are taken from Ref. a. These, together 
with the dileptons from the decays of primary vector mesons, are collectively 
labeled 'background'. It is seen that these background sources describe very 
well the dimuon spectra in the p+W reactions, shown in the figure by solid 
circles. 

However, as can be from the figure, the sum of these background sources 
grossly underestimates the dimuon yield in central S+W collisions, shown in 
the figure by open circles. Since the dimuon spectra are normalized by the mea- 
sured charged particle multiplicity, this underestimation indicates additional 
sources to dilepton production in heavy-ion collisions. This can come from 
the thermalized QGP and/or hadronic phases. So the immediate next step 
is to check whether the contribution from the secondary hadronic processes 
can explain this enhancement. For dilepton spectra at low invariant masses, 
it is well known that the tttt annihilation plays an extremely important role in 
heavy-ion collisions. It is also expected that the other secondary processes will 
play a role in the dilepton spectra in the intermediate mass region. 

Previous thermal rate calculations based on the kinetic theory show that 
in the mass and temperature region relevant for this study, the following sec- 
ondary processes (from the hadronic phase) are very important: ■w' 7r ~ * lh 
up 11, TTLu 11, ttoi -> 11, KK -> 11, and KK* + c.c -> Z/BHtaM Among 
them, the ixa\ — > 11 has been .found to be the most important, mainly because 
of its large cross se^pnE§E!l (Similar conclusion has been drawn for thermal 
photon production EJc3). 

We will concentrate on dilepton production from the secondary processes 
mentioned above. The cross sections for the annihilation of pseudoscalar 
mesons (P) are well known 

87ra 2 fc,„ . „ , , 4m? . . 2m? . , , 

° PP ^m = w i^wi 2 (i - + -^i (3) 

where k is the magnitude of the three-momentum of the pseudoscalar meson in 
the center-of-mass frame, M is the mass of the dilepton pair, and mi is the mass 
of the lepton. It is well known that the electromagnetic form factors \Fp(M)\ 2 
play important role in these processes. The pion electromagentic form factor 
is donimated by the p(770) meson, while that of the kaon is dominated by 
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Figure 4: Left panel: comparison of backgrouds with experimental data in p+W and S+W 
collisions. Middle panel: contributions of various secondary processes to the dimuon spectra 
in central S+W collisions. Right panel: Comparison of the sum of the background and 
secondary contributions to the experimental data in central S+W collisions. 



the phi meson. In addition, at large invarianLpasscs, higher p-like resonances 
such as p(1450) were found to be important Ej. 

The cross sections for up — > 11 and KK* + c.c. ^ 11, which are of the 
pseudoscalar-vector type, have been studied in Ref. cB by fitting to the ex- 
perimental data_for the reverse processes (i.e., hadron production in e + e~ 
annihilation) ql3. High isoscalar vector mesons such as w(1420) and 0(1680) 
play important roles in the electromagnetic form factors of these processes. 
The cross sections from Ref. Ea will be used in this work. The cross section 
for Tiuj — > 11 has the same form as that for 7rp — > 11, but with a different form 
factor, which is_determined from the experimental data for e + e~ — > 7r°7r°7, as 
shown in Ref. C3. 

The cross section for irai — > 11 needs some special attention, since this 
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process has been found toJpe particularly important in the intermediate-mass 
region. Already in Ref. c2l it was shown that the dilepton production rates 
based on different models for the irpai dynamics can_.differ by an order of 
magnitude. This problem was recently revisited in Ref. l3, where a comparative 
study was carried out for both the on-shell properties and dilepton production 
rates using almost all the existing models for the irpa\ dynamics. The results 
were indeed surprising: although most models provide reasonable description of 
the on-shell properties, the corresponding dilepton rate could differ by two-and- 
half orders of magnitude. By using the experimentally-constrained spectral 
function l3, it was found that the effective chiral Lagrangian of Ref. Cj, in 
which the vector mesons are introduced as massive Yang-Mills fields of the 
chiral symmetry, provides the best off-shell, as well as on-shell, properties of 
the irpai dynamics. The cross sections for irai — ► 11 from this model will be 
used in this work. 

The contributions from the secondary processes outlined above are shown 
in the middel Pfpfl of Fig. 4. These are obtained in the relativistic transport 
model of Refs.Jpca, including the HELIOS-3 acceptances, mass resolution, and 
normalization 0. It is seen that the irai process is by far the most important 
source for dimuon yields in this mass region. The ttlo process also plays some 
role in the entire intermediate-mass region. The contributions from tttt, np 
and KK are important around 1 GeV invariant mass. 

In the right panel of Fig. 4, we add the contributions from the secondary 
processes obtained in our transport model to the background, and compare 
again with the HELIOS-3 data for central S+W collisions. It is seen that 
the data can now be nicely explained. Thus we showed for the first time 
the importance of the secondary processes for the intermediate-mass dilepton 
spectra in heavy-ion collisions. This is an important step forward in the use 
of intermediate-mass dilepton spectra as a probe of the phase transition and 
QGP formation. Although the current data do not show any necessity to 
invoke the QGP formation in S-induced reactions, consistent with conclusions 
from the J/^Sf physics, the observation that the secondary processes do play an 
important role in the intermediate-mass dilepton spectra is significant. 

It is appropriate to discuss the issue of the electromagnetic form factor. 
While the form factors for tttt, irp, ttlu, KK and KK* + c.c. are well determined 
from experimental data for e + e~ annihilation, the determination of the 7rai 
electromanetic form factor from the experimental data on e + e~ — > 7r + 7r~7r + 7r~ 
involves some uncertainties. In the previous calculation we assumed that all 
the observed An final state cross section proceeds through the irai intermediate 
state. There are some indication that some of the cross section is dominated 
by the yo(1700) which decays directly to pirn without going through the irai 
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Figure 5: Left panel: comparison of dimuon spectra obtained with the full form factor and 
with the p(700) only. Right panel: comparsion of dimuon spectra obtained with bare and 
in-medium meson masses, respectively. 



intermediate stateEJ. To see the sensitivity of our results on dimuon spectra 
to this uncertainty, we also did a calculation in which the irai form factor 
contains only the normal p(770). The cross section for e + e~ — * irai is then 
very similar to that obtained in Ref. El. The results are shown in the left panel 
of Fig. 5 by the dotted curve. Apparently, with a form factor that excludes 
the p(1700) resonance, the contribution from the nai process is reduced. The 
agreement with the HELIOS-3 data in this case is slightly better. 

Another issue we want to address here is the effects of dropping vector me- 
son masses on the entire dimuon spectra from the threshold to about 2.5 GeV. 
In the previous section we have shown that the enhancement of low-mass dilep- 
tons indicates that the vector meson masses decrease with increase density and 
temperature. This will also affect the dilcpton spectra in the intermediate-mass 
region, through mainly two effects. One is the change of the invariant energy 



12 



spectra of these secondary meson pairs. If the rho meson mass is reduced, then 
the invariant energy of the 7rp collisions should also decrease. The second ef- 
fect enters through the modification of the electromagnetic form factor. Since 
we do know very well how the masses of the p-like resonances change with 
density and temperature, we assume that they experience the same amount of 
the scalar field as the ordinary rho meson, namely, m* , = m PjP i — 2/3g (T (a). 
The results of this calculation are shown in the right panel of Fig. 5. Below 
1.1 GeV and especially from 0.4 to 0.6 GeV, the agreement with the experi- 
mental data is much better when the dropping vector meson mass scenario is 
introduced. At higher masses, the dropping mass scenarios somewhat under- 
estimates the experimental data. In this mass region, however, there might be 
additional contributions from, e.g., secondary Drell-Yan processesE3 that were 
not included in this study. Also, the issues of collisional brodening and baryon 
contribution need to be checked, at these invariant masses. These are under 
study and will be reported on later. 

4 Direct photon: lack of the signal 

Single photon spectra in heavy-ion collisipns at CERN-SPS energies have been 
measured by the WA80 collaboration E3. So far, only the upper bound has 
been determined for the so-called 'thermal' photon spectra, which accounts 
for about 5% of the total observed single photon yield- which is dominated by 
neutral pion decay. In hydrodynamical calculations E3, the absence of signif- 
icant thermal photons has been interpreted as an evidence for the formation 
of a quark gluon plasma. WithouL. phase transition, the initial temperature 
of the hadronic gas found in Ref. t3 is about 400 MeV. This leads to a large 
number of thermal photons from hadronic interaction which is not observed 
experimentally. Including_the phase transition, the initial temperature is low- 
ered to about 200 MeV E3, because of increased degrees of freedom. With 
a lower initial temperature, the thermal photon yield is thus reduced and is 
und to better agree with the WA80 data. However, a more recent analysis 
has shown that if one includes all hadron resonances with masses below 2.5 
GeV the initial temperature can also be lowered to about 200 MeV, and the 
WA80 photon data can then be explained without invoking the formation of 
a quark gluon plasma. This is basically also the conclusion of the detailed 
transport model analysis presented here. 

For the 'thermal' photon spectra we include the decay of p, u>, rj', and a\ 
mesons, as well as two-body processes such as mr — > p"f and ixp — > 717. The 
decay width for p — > inr-f is taken from the model of Ref. E2I, which describes 
well the measured width for p° — > nicy. The u> and a\ radiative decay widths 
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Figure 6: 'Thermal' photon spectra in central S+Au collisions at 200 AGeV from different 
sources. 



are proportional to Ip^l, with the coefficients determined from the measured 
width. For the two-body cross sections, we use the results ofE3, which do 
not include the contribution from an intermediate a\ meson. The latter has 
already been included in our model as a two-step process. 

The contributions to the so-called 'thermal' photons from the decay of 
mesons and baryons, and from two-body scattering are shown in Fig. 6. The 
u) radiative decay is found to be the most important source for photons with 
transverse momenta above about 0.5 GeV. The contribution from the a\ ra- 
diative decay is somewhat larger than_that from direct -np — > ttj. This is in 
agreement with the conclusion of Ref. c3 that including a\ the 7rp contribution 
increases by about a factor 2-3 in the relevant temperature and photon energy 
region. 

The contributions from rf and S° radiative decays are restricted to photons 
with transverse momenta below 0.5 GeV, because the mass differences between 
hadrons in the intial and final states in these processes are small. Photons with 
transverse momenta below 0.2 GeV come chiefly from the decay of S° and tttt 
scattering. The reaction tttt — > is endothermic, with most of the available 
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energy going into the rho meson mass. This cross section actually diverges as 
the photon energy goes to zero. We have included a low-energy cut-off as in 
Ref. E3. The choice of this cut-off parameter affects basically only the photons 
with transverse momenta below about 0.1 GeV, where no experimental data 
are available. 
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Figure 7: 'Comparison of 'Thermal' photon spectra in central S+Au collisions at 200 AGeV 
from different calculations. 

The importance of the radiative decay contribution from primary omega, 
mesons was not seen in thermal rate calculations such as those of Ref. E3. 
Non-equilibrium effects and dynamical evolution of the collision attribute to 
the difference between the thermal and transport model calculations. As the 
colliding system expands, the two-body contribution becomes less important 
as compared with that of one-hody, since the former is proportional to density 
squared. It was found in Ref. E3 that in the first a few frn/c, when the system 
is the densest, the contribution from np scattering is indeed more important 
than that from u> decay, in agreement with the findings of the thermal rate 
calculation in Ref. EE However, the rate from irp damps with time very fast. 
After 3-5 frn/c, the emission rate from uj meson becomes more important. 
Adding on top of this the contribution from the omega meson decay at the 
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freeze out, it is easy to understand why the u) meson radiative decay becomes 
as important as that from irp collisions in a dynamical calculation. 

The WA80 single photon data have been looked at by several groups, 
mostly based onJ^drodynamicaljnodcls In Fig. 7, we show our results together 
with those fromE3 andcJ. In Ref. Ej several different equations of state with and 
without quark gluon plasma formation were considered. Shown in the figure 
by dotted line is their results based on EOS H. This equation of state included 
about the same number of hadron degrees of freedom as in our transport 
model. It is very interesting to see that their results are very similar to ours, 
although the dynamical models used are quite different. In Ref. c3 two different 
equations of state, one with quark gluon plasma formation and one with pure 
pionic gas, were considered. Shown in the figure is their results based on 
the pionic gas equation of state. Because of very limited number of hadronic 
degrees of freedom in this equation of state, they needed a very high initial 
temperature to account for the final observed hadron abundances. This led to 
large photon yield. 

5 SUMMARY AND OUTLOOK 

In summary, we have studied dilepton and photon production from both proton- 
nucleus and nucleus- nucleus collisions using the relativistic transport model 
with initial conditions determined by string fragmentation from the initial 
stage of the RQMD model. It is found that the dilepton spectra in proton- 
nucleus reactions measured by the CERES and the HELIOS-3 collaboration 
can be well understood in terms of conventional mechanisms of Dalitz decay, 
direct vector meson decay, decay of charmed mesons, and initial Drell-Yan pro- 
cesses. For dilepton spectra in central heavy-ion collisions, these conventional 
mechanisms, however, fail to explain the data, in the entire mass region from 
threshold to about 2.5 GeV. 

Including the contribution from pion-pion annihilation, which is important 
in the mass region from 2m 7r to m PjW , removes some of the discrepancy. But 
the theoretical prediction is still substantially below the data in the low mass 
region and somewhat above the data around m pM - The theoretical results are 
brought into good agreement with the data when reduced in-medium vector 
meson masses are taken into account. In the intermediate-mass region, we have 
included additional secondary processes involving meson-meson collisions. The 
cross sections for these processes are constrained by the experimental data for 
e + e~ annihilation. We found that these secondary processes, especially the 
7rai — > 11, play very important role in the intermediate-mass region. We have 
also calculated, within the same dynamical model, the thermal photon spectra 



1G 



in central S+Au collisions. In both free meson masses and in- medium meson 
masses scenarios, our results do not exceed the upper limit deduced from the 
experiments by the WA80 collaboration. 

So far, we have not included explicitly collisional brodening in our calcu- 
lation. Since the rho meson is treated as a dynamical particle in our transport 
model, its collisions with meson (mainly pions) and baryons (mainly nucleons) 
are included. This might reflect, at least partly, its collisional brodening in 
hot and dense matter. On the other higher resonances, such as p(1450) and 
p(1700), are not treated dynamically. Their effects are included through elec- 
tromagnetic form factors. Therefore, collisional brodening is not considered for 
these higher resonances. Since they usually have a natural width of 200-300 
eV, a collsional brodening width of about 50 MeV for these higher resonances 
shall not affect our results dramatically. This, however, needs to be dealt 
with quantitatively. 

The current investigation can be extended to higher incident energies, such 
as those of RHIC collider, by combining the cross sections (or thermal rates) 
obtained in this study with, e.g., hydrodynamical models for the evolution of 
heavy-ion collisions at the RHIC energies. This kind study will be very useful 
for the determination of hadronic background in the dilepton spectra, and for 
the clear identification of dilepton spetra from the QGP. 
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